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Planar tunneling junctions of the presumably unconventional superconductor FeSe were prepared 
and investigated. The junctions consisted of FeSe/AlOji/Ag multilayers patterned lithographically 
into mesa structures. Bias voltage dependent differential conductivities dI/dV(V) of junctions in the 
tunneling, intermediate barrier strength, and metallic regimes were investigated. Depending on the 
barrier type between two and four features of the conductivity were obtained, which are discussed in 
the framework of multiple superconducting energy gaps. Specifically we reproduced with all barrier 
types an established energy gap and discovered an additional large gap structure. From two further 
conductivity features presumable one originates from another superconducting gap and the other 
from resonant states. 
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I. INTRODUCTION 

Immediately after the discovery of the first iron based 
superconductor [I| , many publications discussed possible 
unconventional pairing mechanisms [H- These mech¬ 
anisms are typically associated with specific symmetries 
of the order parameter. An extended s-wave state was 
proposed by Mazin et al. Q, which is supported by re¬ 
cent theoretical work @ . Additionally, also d-wave states 
were predicted Q- However, only scarce direct experi¬ 
mental evidence for such symmetries is not available up 
to now. 

Amongst the iron based superconductors FeSe has the 
simplest crystal structure with a critical temperature of 
Tc — 8.5 K 0. Critical temperatures up to 11.4K were 
obtained by biaxial strain resulting in an elongation of 
the c-axis [l0|| while a compression of this axis suppresses 
superconductivity [ll|. The critical temperature can be 
further increased up to Tc — 37 K by applying hydro¬ 
static pressure of 4.1 GPa Q. 

Theoretical descriptions of the electronic properties of 
FeSe are usually based on a five-band model with elec¬ 
tron pockets at (0, ±7r) and (±7r, 0) and hole pockets at 
(0,0) [l^. On this multiband Fermi surface different 
types of extended s-wave and d-wave order parameters 
are discussed, which include nodal as well as nodeless 
superconducting energy gaps. In-situ scanning tunnel¬ 
ing spectroscopy (STS) on MBF grown FeSe thin films 
was interpreted as evidence for gap nodes by Song et 
al. [l^. These authors observed one very clear single V- 
shaped energy gap with As ss 2.2meV [T3l - [l^ . which 
corresponds to ratios of 2AB/kBTc = 5.5 — 6.4. An ad¬ 
ditional conductivity feature at 2Ac ~ QAkBTc, was in¬ 
terpreted as strong coupling to a spin fluctuation mode 
(l7l| . However, based on the multiple Fermi surfaces of 
FeSe, also multiple superconducting energy gaps are pos- 
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sible. Recent STS investigations of high quality FeSe 
single crystals by Kasahara et al. observed addition¬ 
ally to a clear V-shaped energy gap with Ab = 2.5 meV 
{2Ab/ kBTc = 6.8) a shoulder feature in the differential 
conductivity at ~= 3.5 meV, which they interpreted as a 
second superconducting gap with 2Ac/kBTc = 9.5 [T8l| . 
Also by Andreev spectroscopy on break-junctions mul¬ 
tiple features were observed and explained by two en¬ 
ergy gaps with A^ w 1.0 meV {2AaJ_^Tc = 1.9) and 
Ab ^ 2.6 meV {2AB/kBTc = 6.7) [l^. Summarizing 
the available literature about the superconducting en¬ 
ergy gap of FeSe agreement about a clear gap feature 
with 2 Ab//cbTc = 5.5 — 6.7 exists. However, the ob¬ 
servation of additional conductivity features and their 
interpretation is inconsistent. 

Here we present spectroscopic measurements of the 
bias voltage dependent differential conductivity of pla¬ 
nar FeSe/A10a;/Ag tunneling junctions. We discuss three 
different kinds of spectra including tunneling like behav¬ 
ior, metallic point contact like behavior and data sets 
showing characteristics of both junction types. 


II. JUNCTION PREPARATION AND 
MEASUREMENT PROCEDURE 

Epitaxial FeSe thin films were prepared by rf-sputter 
deposition from two elementary 2 inch targets (Fe: 
Mateck GmbH, purity 99.99%; Se: Lesker, purity 
99.999%) using AJA International A320-XP-MM sput¬ 
tering sources. Superconducting thin films in (001)- 
orientation were obtained on MgO(lOO) substrates. For 
more information on the growth process and structural 
and morphological characterization of the /3-FeSe thin 
films please see [^ . 

For junction preparation, on top of the /3-FeSe thin 
films (thickness d Ri 50 nm) thin Al layers (2.8 nm) were 
sputter deposited and plasma oxidized for tunneling bar¬ 
rier formation. Sputter deposited Ag thin films (d = 
50 nm) served as counter electrodes (d = 50 nm). Finally, 
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the three layer samples were patterned into mesa struc¬ 
tures of different sizes (from 60x60p,m^ to lOOxlOO/rm^) 
by optical lithography and ion beam etching. 

The bias voltage dependent differential conductiv¬ 
ity of the junctions was measured by a standard a.c.- 
modulation technique (Vac = 50 /iV) in a ^He-cryostat. 


III. SPECTROSCOPIC RESULTS 

None of our junctions showed a differential conductiv¬ 
ity dI/dV{V) approaching zero at zero bias, indicating 
that either the FeSe thin films were never homogeneously 
superconducting across the junction area or the AlOa, 
thin layer never represented a perfect tunneling barrier. 
Figure 1 shows the differential conductivity of a junction, 
which is closest to tunneling dominated behavior as in¬ 
dicated by the largest ratio of the conductivities outside 
and inside the gap. 



FIG. 1: Differential conductivity dI/dV„orm{V) of a 

FeSe/A10a;/Ag junction at T — 0.3 K, 2K, 4K, and 5K, nor¬ 
malized to dI/dV{V) measured at T = 6K. The inset shows 
the conductivity at T = 0.3 K and 6K without normalization. 


The inset of Figure 1 shows the junction conductiv¬ 
ities obtained at T = 0.3 K and at the temperature 
T = 6 K {Tc ~ 5 K) to which the conductivity curves dis¬ 
played in the main panel were normalized. The normal¬ 
ized spectrum measured at T = 0.3 K shows a V-shaped 
gap feature in an energy range of ±1.2 meV around zero 
bias. Additionally, a second gap feature with conductiv¬ 
ity peaks at ±3.7 meV is obvious. Identifying both fea¬ 
tures with superconducting energy gaps 2AB/kBTc = 5.6 
and 2AB/kBTc = 17.2 are obtained. 

A tunneling barrier with pinholes can result in a 
large variety of structures in the differential conductiv¬ 
ity. These structures are often caused by local heat¬ 
ing effects. However, small pinholes can behave anal¬ 
ogous to point contacts representing metallic contacts 



FIG. 2: Differential conductivity dl/dVnorm(V) of a 
FeSe/A10a;/Ag junction at T = 0.3 K, normalized to 

dI/dV(V) measured at T = 6K and BTK-fit assuming two 
energy gaps (Z = 0.3). The right inset shows the conductiv¬ 
ities at T = 0.3 K, 3K and 5K without normalization. The 
left inset shows the differential conductivities at T = 0.3 K, 
2K, and 3K. 


with ballistic electron injections resulting in Andreev re¬ 
flection as described by the BTK-theory [2lj. In Fig¬ 
ure 2 a normalized Andreev-like spectrum obtained mea¬ 
suring a FeSe/A10a;/Ag junction is presented for tem¬ 
peratures between 0.3 K and 5K. The critical temper¬ 
ature of this sample amounts to Tc = 5 ± 0.5 K, as 
indicated by an almost constant dI/dV{V) curve mea¬ 
sured at this temperature. As shown in the right inset 
of Figure 2, the dI/dV{V) curve obtained just above 
Tc is shifted to lower conductivities compared to the 
curves measured below Tc. This indicates that in the 
normal state of the FeSe base electrode, the pinhole(s) 
and the FeSe bulk resistivity within the mesa structure 
both contribute about half to the total mesa conductiv¬ 
ity. Similar effects were observed in other ironpnictide 
superconductor-normal metal-superconductor junctions 
[22j| . Another deviation of the junction properties from 
the conventionally expected behavior of ballistic point 
contacts is the almost missing temperature dependence 
of the experimental differential conductivity up to Tc. 
The very weak temperature dependence observed is dis¬ 
played in the left inset of figure 2. However, following 
the procedure of [1^, the experimental dl/dViV) curve 
measured at 0.3 K can be fitted very well using the BTK- 
theory with Dynes broadening as shown in the main panel 
of Figure 2. The magnitudes of the gaps used in the 
fit are As = 1.2 ± 0.2 meV and = 4.5 ± 0.2 meV 
corresponding to ratios of 2AB/kBTc = 5.6 ± 0.9 and 
2Az,/fcBTc = 20.9 ±0.9. 

The third type of spectra obtained investigating 
FeSe/A10a;/Ag junctions, shows characteristics resem¬ 
bling both tunneling and Andreev reflection, which is 
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typical for junctions with intermediate barrier strength 
(2l|. Compared to the junction with the Andreev-like 
spectrum discussed above, only a small contribution 
5%) of the bulk resistivity of the FeSe to the total 
measured mesa structure conductivity is deduced from 
the shift to lower conductivities in the normal state. In 
Figure 3 the temperature dependent conductivity of a 
junction with several features, which could be associated 
with superconducting energy gaps, is shown. 



FIG. 3: Differential conductivity dI/dV„orm{V) of a 

FeSe/A10a:/Ag junction at T = 0.3 K, 1K,2K, ... 8K, nor¬ 
malized to dI/dV{V) measured at T = 9K. The inset shows 
the conductivity at T = 1.3 K without normalization, labeling 
different dI/dV(V) features. 


The broad hump feature B at ±1.9 meV can be asso¬ 
ciated with the energy gap Ab discussed above (Figs. I 
and 2), corresponding to 2AB/kBTc = 5.6 ± 1.2. The 
broad energy gap Ap discussed above corresponds to 
a kink labeled D in the dI/dV(V) curve at ±7.5 meV 
{2Ab)/ kBTc = 21.5 ± 2.9). Additionally, there are 
four more features of the differential conductivity, from 
which the very sharp feature A appears at ±1.1 meV 
{2TA/kBTc = 3.2 ± 0.6). At higher bias voltages a clear 
shoulder at rs 3.7 meV (labeled C) is observed, corre¬ 
sponding to 2Ac/kBTc = 10.6 ± 1.2 meV. Further fea¬ 
tures are visible at ±0.4 meV and zero bias, which, based 
on the study of Fe atoms on a FeSe surface by Song et 
al. [l^, can be explained as resonance states. 

Please note that the features discussed above are repro¬ 
duced by several junctions, which are in the intermediate 
barrier strength regime. The differential conductivity of 
another junction in the same regime is shown in Figure 
5. Although the shape of the features discussed above is 
partially different and the critical temperature is strongly 
reduced to 4 K, all potential gaps can be found again at 
the same 2A/kBTc values. 

The features A-D appear at bias voltages of 0.3 ± 
O.lmeV, 0.9±0.1meV, 3.4±0.4meV, and 5.5±1.0meV, 
corresponding to 2TA/kBTc = 1.9 ± 0.6, 2AB/kBTc = 



FIG. 4: Differential conductivity dI/dVnorm(V) of a 

FeSe/A10a;/Ag tunneling junction at T = 0.3K, 2K, 3K, and 
4K, normalized to dI/dV{V) measured at T = 5K. The inset 
shows the conductivity at T = 1.4K and 5K without normal¬ 
ization. 


5.4 ±0.6, 2Ac/kBTc = 9.9 ± 2.3, and 2AD/fcsFc = 
16.0 ±5.8. 


IV. DISCUSSION 

In the differential conductivity of all junctions two to 
four features can be associated with superconducting en¬ 
ergy gaps. As such, they should scale with the criti¬ 
cal temperature, which allows the comparison of sam¬ 
ples with different quality reflected by different Tc val¬ 
ues. Thus in Table 1 an overview of the 2A/kBTc values 
is given and compared with previously published data. 

It is obvious that in all data sets an energy gap with 
2AB/kBTc = 5.4 —6.8 is observed. Additionally, also the 
shape of the dI/dV(V) curves confirm that feature B is 
a clear superconducting energy gap. 

Concerning the other gap like features none of them is 
observed simultaneously in all data sets. We observed a 
sharp conductivity feature at 2T A/kBTc = 1.3 — 3.8, an 
energy range in which the existence of a superconducting 
energy gap was proposed based on break-junction spec¬ 
troscopy by Ponomarev et al. [H. However, the features 
A in the differential conductivity spectra of our junc¬ 
tions with intermediate barrier strength are very sharp 
and more strongly suppressed by increasing temperature 
compared to the features B. In this sense they resemble 
resonance states as already discussed above for the lowest 
energy and zero bias conductivity peaks shown in Fig. 3. 

Our data sets support the existence of a conductivity 
feature at 2Ac/kBTc = 9.2 — 11.8, which was proposed 
by Kasaharaet al. based on STS to originate from a su¬ 
perconducting energy gap [l^. However, considering the 
shape of the dI/dV(V) curves (Figs. 3 and 5) this could 





















4 


from Figure 

2T A/kBTc 

2Ab/ kBTc 

2Ac/ksTc 

<1 

1 

- 

5.6 ±0.9 

- 

17.2 ± 1.9 

2 

- 

5.6 ±0.9 

- 

20.9 ±0.9 

3 

3.2 ±0.6 

5.5 ± 1.2 

10.6 ± 1.2 

21.5 ±2.9 

5 

1.9 ±0.6 

5.4 ±0.6 

10.4 ± 1.2 

19.7 ±2.3 

Kasahara et al. 

- 

6,8 

9,5 

- 

Song et al. [W, ^ 

- 

5, 5-6, 4 

9,4 

- 

Ponomarev et al. [1^ 

1,9 

6,7 

- 

- 


TABLE I: Overview of 2A./kBTc for all potential energy gaps / conductivity features. 


also be strong coupling feature as proposed by Song et 

ai.[i3. 

Most important, we measured in all of our junctions 
an additional gap feature at relatively high energies with 
2IS.D/ksTc = 15.3 — 24.3, which cannot be found in the 
existing STS or break-junction spectroscopy data of FeSe. 
However, this feature appeared clearly in the most tun¬ 
neling like (Fig. 1) and in the most metallic (Fig. 2) junc¬ 
tion and shows in both cases the characteristics expected 
for a superconducting gap. 

As a simple test of the identification of the conductivity 
features B and C with superconducting energy gaps their 
temperature dependence is plotted in Figure 5. 



T/Te 


FIG. 5: Temperature dependence of the normalized energy 
gaps, As and Ac as well as of the feature Ta from Figure 3. 
Additionally, the temperature dependence Ascs(T')/A(0) as 
given by the BCS theory is plotted. The inset shows Ab{T), 
Ac{T), and rA(T) without normalization. 


V. SUMMARY 

The superconducting state of FeSe was investigated by 
tunneling spectroscopy on planar FeSe/A10a;/Ag junc¬ 
tions. By an evaluation of the bias voltage dependent 
differential conductivities of junctions in the tunneling, 
intermediate barrier strength, and metallic regimes 
common features could be identified and associated with 
the superconducting density of states. Specifically, four 
characteristic features were observed. Whereas two of 
them show clear properties of superconducting gaps with 
2AB/kBTc = 5.4 - 6.8 and 2AD/kBT^ = 15.3 - 24.3, 
the interpretation is less clear for the other two. There 
are indications of a strong coupling effect with a boson 
excitation with uJboson ~ l.SmeV as previously pro¬ 
posed by Song et al. [II- The fourth feature with 
2VA/kBTc = 1.3 — 3.8 is presumably an impurity in¬ 
duced resonance state in the superconducting energy gap. 

Financial support by the DFG (Jo404/6-l) is acknowl¬ 
edged. 


As well Ab{T) as Ac{T), for which the association 
with an energy gap or a strong coupling effect is ambigu¬ 
ous, follow the prediction of the BCS theory supporting 
their identification as energy gaps. Additionally, also the 
temperature dependence of the resonance feature A is 
plotted, which disappears already well below Tc- 
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